Completeness of the classical 2D Ising model and universal quantum computation 
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We prove that the 2D Ising model is complete in the sense that the partition function of any 
classical g-state spin model (on an arbitrary graph) can be expressed as a special instance of the 
partition function of a 2D Ising model with complex inhomogeneous couplings and external fields. In 
the case where the original model is an Ising or Potts-type model, we find that the corresponding 2D 
square lattice requires only polynomially more spins w.r.t the original one, and we give a constructive 
method to map such models to the 2D Ising model. For more general models the overhead in system 
size may be exponential. The results are established by connecting classical spin models with 
measurement-based quantum computation and invoking the universality of the 2D cluster states. 

PACS numbers: 03.67.-a, 03.67.Lx, 75.10.Hk, 75.10.Pq, 02.70.-c 



1. — Introduction. Classical spin models such as the 
Ising and Potts models are widely studied in statistical 
physics, as they provide important toy models for mag- 
netism and as they can be mapped to numerous inter- 
esting problems in physics and mathematics [3, [2] . The 
geometry of a model, in particular its spatial dimension, 
plays an important role with respect to the physical prop- 
erties of the system and the possibility of finding (approx- 
imate) solutions. For instance, it is known that evalua- 
tion of the partition function of the Ising model with 
magnetic fields is easy in ID, while on a 2D square lat- 
tice this problem is already NP-hard Q . 

In this paper we study the interrelations between 
classical g-state spin models on different geometries (or 
graphs), and find that the 2D Ising model (which has 
q = 2) plays a distinguished role in this study. We con- 
sider mappings that leave the partition function — and 
hence all thermodynamical quantities, such as free en- 
ergy or magnetization, derived from it — invariant (see 
also 0, HI). As the main result of this paper, we prove 
that the 2D Ising model is complete in the sense that 
the partition function of any classical g-state spin model 
can be expressed as a special instance of the partition 
function of a 2D Ising model with inhomogeneous cou- 
plings. More precisely, given a partition function Zq of 
a g-state spin model on an arbitrary graph — which may 
be, e.g., a lattice of arbitrary dimension or involve long- 
range interactions — there exists a 2D square lattice of 
enlarged size, and suitably tuned nearest-neighbor cou- 
pling strengths and magnetic fields, such that the parti- 
tion function of the Ising model on this lattice specializes 
to Zq. Furthermore, in the case where the original model 
on the graph G is an Ising or Potts-type model, we find 
that the corresponding 2D square lattice requires only 
polynomially more spins w.r.t the original one. For more 
general models the overhead in system size may be ex- 
ponential. However, one important remark needs to be 
made: in order to achieve this result, one has to allow 
for complex couplings in the 2D partition function — thus 



leaving the "physical" regime of the model. 

The results are proven by relating the problem at 
hand to insights from quantum information theory, more 
particularly to the area of measurement-based quantum 
computation (MQC). The latter is a recently established 
paradigm for quantum computation where quantum in- 
formation is processed by performing sequences of single- 
qubit measurements on a highly entangled resource state 
Q. In order to obtain our results, we first prove that 
the Ising partition function on an arbitrary graph (with 
external field) can be written as the overlap between an 
entangled quantum state and a complete product state — 
thus generalizing a construction which we introduced in 
Ref. [Bj]; see also Ref. Q- This formulation allows us 
to make a connection with MQC. In particular, we prove 
that the entangled state corresponding to the Ising model 
on a 2D square lattice, is (a variant of) the 2D cluster 
state Q . The latter is known to be a universal resource 
state for MQC in the sense that every quantum state 
can be obtained by performing a suitable sequence of 
single-qubit measurements on a sufficiently large 2D clus- 
ter state. This quantum universality feature of the 2D 
cluster states leads to the result that the 2D Ising model 
is complete in the sense specified above. 

2. Classical Ising model. — We consider the classical 
Ising model involving N two-state spins (si, S2 ■ ■ . sjv) = 
s, where s a = ±1. The spins interact pairwise ac- 
cording to an interaction pattern specified by a graph 
G = (V, E) with vertex set V and edge set E, and 
the coupling strengths are denoted by J a (,. Moreover, 
the spins are subjected to local magnetic field terms h a . 
The Hamiltonian of the system is given by Hq(s) := 



b}£E 



JabSaSb 



In other words, we 



consider a general inhomogeneous Ising model on an ar- 
bitrary graph. The partition function Zq is defined by 
Z G {{J abl h a }) := Y. e~^ s \ where = (k B T)~\ with 
ks the Bolzmann constant and T the temperature. 

3. Quantum formulation. — We now show how the 
partition function Zq can be expressed in a quantum 
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physics language. Let G be the graph with n = \V\ + 
\E\ vertices and 2\E\ edges which is obtained from G by 
adding at each edge {a, 6} £ E an additional vertex ab 
and thus "splitting every edge in half" (see Fig. [1]). We 
will call G the decorated version of G. The vertex set of G 
is thus given by the union of the original vertex set V of G 
and the set Ve — {ab | {a, b} £ E} corresponding to edges 
of G — note that we label the vertices in Ve by double 
indices, indicating their origin in the corresponding edge 
of G. We now consider an n-qubit state \<Pq), defined on 
a set of qubits labelled by V U Ve, which is defined to 
be the graph state 0, associated with the decorated 
graph G. In particular, \iPq) is the (unique) joint fixed 
point of the \V\ + \E\ stabilizing operators K a and K a b, 



K a = I*"' JJ X iab ^ 

b:{a,b}£E 

K ab = Z {ab) Z {a) Z {b \ 



(1) 



for every a G V , and for every e = {a, b} £ E. Here 
X and Z denote the Pauli spin matrices, and the upper 
indices indicate on which qubit is acted. 

We can now express the partition function as follows: 



Z G ({ J ab ,h a }) = 2\ y \' 2 -{a\<pQ). 
In this expression, 



\a) = 



\ab£V E 



\ot a b) \a a ) 



(2) 



(3) 



is a complete product state specifying the coupling 
strengths of the Ising model. In particular, \a a b) = 



|0> 



-0Ja 



|1) is an (unnormalized) one-qubit state 



(acting on qubit ab) determined by the interaction 
strength between particles a and b. Similarly, |a a ) = 



|0> 



1} is an (unnormalized) one-qubit state 



(acting on qubit a) determined by the local magnetic 
field at particle a. Expression |(3J) shows that Zq can be 
obtained by calculating the inner product of the graph 
state \<Pq) and a complete product state. The choice of 
the product state allows one to specify the couplings of 
the Hamiltonian and the temperature, while the struc- 
ture of the graph state reflects the interaction pattern. 

To show that Eq. @ holds, we use that \<p@) can be 
written as \^>q) oc ^ t \B T t)\t), where t is a binary vector 
of length |V|, B is the incidence matrix of the graph G, 
and by writing out the sum (a|^g). The construction 
of \v>q) can be viewed as a generalization of the one we 
introduced in Ref. While in Ref. [B[ each qubit 

was associated with an edge of the graph G, here we 
have two types of vertices: one subset Ve associated to 
edges ( "edge-qubits" ) and one to vertices V ("vertex- 
qubits" ) . This enlarging of the system size allows one to 
treat also local terms in the Hamiltonian (whereas Ref. 
[B[ only dealt with zero external field). In addition, the 



stabilizer of the state \<Pq) can be immediately obtained 
from the graph G describing the interaction pattern (or 
its decorated version G), as in Eq. |T]) and Fig. Q] 



G 



G 



FIG. 1: (Color online) Decorated graph G (right) correspond- 
ing to a 2D lattice G (left). Green (dark) dots indicate ver- 
tices originating from the vertices V of G, while red (light) 
dots indicate vertices Ve originating from edges E of G. 



4- MQC and the 2D cluster states. — We now turn our 
attention to measurement-based (or: "one-way") quan- 
tum computation, and establish a relation to the parti- 
tion function of the 2D classical Isingmodel via Eq. @. 

The one-way quantum computer [6| is a recently de- 
veloped model for quantum computation, where compu- 
tations are realized by performing single-qubit measure- 
ments on a highly entangled substrate state called the 
2D cluster state \C) [!}; the latter is a graph state (loj 
associated to a 2D square lattice C. 

A particular feature of the one-way quantum computer 
is that it is universal. This means that any n-qubit quan- 
tum state can be prepared, up to local unitary Pauli op- 
erations, by performing sequences of single-qubit mea- 
surements on a d x d cluster state \C) of sufficiently large 
system size M = d 2 . This property of the 2D cluster 
states immediately implies that every n-qubit quantum 
state \tp) can be written in the following way: 



= 2( M -™)/ 2 (7®(/3|)|C). 



(4) 



This formula represents one "measurement branch" of a 
one-way computation performed on an M-qubit cluster 
state, yielding the state \tp) (up to a local operation S) 
as an output state on the subset of qubits which has not 
been measured. The dual product state (f3\ = J -(/9_ 7 '|, 
which acts only on the measured qubits, is determined 
by the bases and the outcomes of the different steps in 
the computation. The local unitary operator £ ("cor- 
rection operator") acts on the unmeasured qubits (i.e., 
on the Hilbert space of the tensor factors of S are 
always instances of Pauli operators: Si e {I, X, Y, Z}, 
The prefactor 2( M- ™)/ 2 reflects the fact that the success 
probability of every measurement branch is 2 n ~ M . 

As proved in Ref. 0], for all n-qubit states that 
can be efficiently prepared in the circuit model, i.e., by a 
polynomial sequence of two-qubit gates, the required size 
M of the cluster state in Eq. (j4|) scales polynomially with 
the number of qubits: M oc poly(n). Moreover, in this 
case the measurement bases \f3j) as well as the correction 
operations S can be efficiently determined. Since any 
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graph state on n qubits can be prepared using at most 
0(n 2 ) controlled-phase gates (lo| . it follows that an arbi- 
trary n-qubit graph state can be written in the form 
P|) with M — poly (n). Furthermore, for the preparation 
of graph states every single-qubit state \/3j) can always 
be chosen to be one of the X-, Y- and Z-eigenstates. 

Also l^j) (i.e., the state \<p G ) where G = C is the 2D 
square lattice) is a universal resource. This is because the 
2D-cluster state \C) can be deterministically generated 
from \<fic) (up to a local correction) by performing single- 
qubit Y-measurements on all qubits in Ve- This fact was 
already noted in Q. As a consequence, one has 

E'|C) = 2l £ l/ 2 (/®(0 F r-)|^), (5) 

where \0y) Ve is a tensor product of the (+l)-eigenstate 
of Y on all edge-qubits, and £' is a local correction. 

5. Universality of the 2D Ising model. — We are now 
ready to establish the connection between the evaluation 
of Ising partition functions and universal MQC. To this 
aim, consider the Ising model on a graph G. The parti- 
tion function Z G can be expressed in the form @. Now 
consider the following procedure. 

First, the graph state \<p G ) is written in the form (|4|) 
when taking \ip) = \<p G ). Together with Eq. (|5|), this 
implies that the partition function Z G can be written as 

Z G ({J ab ,h a })=A-{ 7 \<p d ), (6) 

where A is a constant and I7) is a product state, I7) = 
T,\a)(g>^'\P)(g>\0 Y ) VE ■ Note that, as \ip G ) is a graph state, 
the system size of the 2D cluster state grows polynomially 
with the size of G. Furthermore, |/3) consists of X-, Y- 
and Z-eigenstates. 

Now, applying Eq. ([2]) to the 2D Ising model, the 
overlap between \<pg) and a complete product state cor- 
responds to a 2D Ising partition function Z 2 £>, evaluated 
in certain couplings {J[y h' { } determined by I7). This 
allows us to conclude that Zq can be written as follows; 

Z G ({J ab ,h a }) <x Z 2D ({Jl j ,h' i }). (7) 

In other words, the Ising partition function on an arbi- 
trary graph can be recovered as a special instance of the 
Ising partition function on a 2D square lattice. 

Note that, in the above sequence of arguments, one 
step is particularly crucial, namely the universality of 
the 2D cluster states: this property is used to "map" an 
arbitrary state \<p G ), and hence the associated partition 
function, to the 2D cluster state, i.e., all states can be 
"reduced" to this single structure. 

We give a few remarks regarding this construction. In 
Eq. ([6]), note that the product state I7) is determined by 
both the interaction graph G and the couplings { J a b, h a } 
of the original model. For, on the one hand, it contains 
the states \a a b) and \a a ) encoding the couplings of the 
original model; on the other hand, I7) contains states |/3y) 



and |0y) corresponding to the sequence of one-qubit mea- 
surements which are to be implemented in order to gener- 
ate \<p G ) from the universal resource |<£<j). In going from 
Eq. ([6]) to Eq. ([7|), the state I7) in turn determines the 
couplings in which the 2D model is to be evaluated. Note 
that the decorated cluster state \<pg) has vertex-qubits 
and edge-qubits. The factors of I7) acting on the edge- 
qubits determine the pairwise interactions J 8 '-, whereas 
the factors of I7) acting on the vertex-qubits determine 
the external fields h[. The tensor factors of (7) which 
act on the edge-qubits are all equal to |0y) cx |0) + i\l). 
This implies in particular that, in ((7|), only homogeneous 
pairwise couplings Jy need to be considered. Further- 
more, due to the imaginary unit "i" in |0y), these cou- 
plings generally lie in a complex parameter regime; in 
particular, one can show that /?J|- = —iir/4 is a cor- 
rect choice. Also, the fact that the can be chosen to 
be homogeneous implies that all information regarding 
the pairwise couplings 3 ab and external fields h a of the 
original model, and the graph G of this model, will be 
encoded in the factors of I7) acting on the vertex-qubits, 
and thus in the external fields h\ (which will typically 
be inhomogeneous). We further remark that the part 
of 1 7) acting on the vertex-qubits generically also cor- 
responds to complex interaction strengths h\ (e.g., \(3) 
may contain Y-eigenstates) . A special role is played 
by those factors of \(3) which are equal to Z-eigenstate 
0) oc e oc |0) + e _oo |l). These states give rise to "infinitely 
large" external fields at the corresponding vertices, which 
effectively corresponds to a boundary condition. 

In conclusion, the universality of the 2D cluster states 
\(pg) in the context of MQC, implies that the Ising par- 
tition function on any graph can be expressed as a 
special instance of a (polynomially enlarged) 2D Ising 
model with complex, homogenous pairwise interactions 
and complex, inhomogenous external fields. Note that 
even though such complex interaction strengths do not 
correspond to physical models, considering the partition 
function as a function with complex arguments is com- 
monly done, e.g., in the context of evaluating the Tutte 
polynomial or finding (complex) zeros of Z G to identify 
phase transition points 0]. 

6. Generalizations to q-state models. — Our results 
can also be generalized to g-state spin models such as 
the Potts model We showed in Ref. [B| that the 
partition function of a g-state Potts model on a graph 
G = (V, E) can be written as the overlap between a 
stabilizer state \tp%) and a complete product state \x). 
Z G oc {x\^P%)- Similar to the treatment of the Ising 
model, the state |</r=,) depends only on the graph, and 
the state \x) = ® Qb \Xab) is a complete product state 
depending only on the couplings of the model. However, 
the main difference is that the single-particle systems are 
no longer qubits, but g-dimensional systems. E.g., one 
finds \xab) = e /3Ja6 |0) + J2k=i 1^) B' Interestingly, the 
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partition function of such a g-state model (for arbitrary 
graphs) can again be expressed as a special instance of 
the partition function of the 2D-Ising model (with q = 2) 
and complex parameters — again using the connection to 
MQC. To achieve this this, we use that any g-dimensional 
product state can be mapped by a suitable unitary op- 
eration to a product state of m q — [log 2 q\ qubits; e.g., 
\Xab) = Ul b \0)® m «. As q is fixed, the unitary U a b can 
be prepared with a constant number of two-qubit gates. 
Being a stabilizer state, \<p^) is preparable by a poly- 
sized (qubit) circuit. It follows that Zq can be written 
as the inner product of an efficiently preparable state 
\<p) := < S> ab U a b\tp q ^,) (which is now regarded as a multi- 
qubit state) with a product state \0)® m i\ E \. The uni- 
versality of \tpg) for MQC now implies that \<p) can be 
obtained by performing single-qubit measurements on a 
polynomially enlarged cluster state \<fc). In particular, 
Eq. (j4|) can be applied to = \ip). Using a similar 
argument to Section 5, this implies that the Potts model 
partition function is a special instance of the partition 
function of a polynomially enlarged 2D-Ising model with 
properly tuned complex parameters and two-state spins. 

The above strategy can even be applied to g-state mod- 
els beyond the Potts model, e.g. to all models on directed 
graphs where the Hamiltonians are arbitrary functions 
of the difference (modulo q) between spin values, includ- 
ing arbitrary local terms, while still obtaining a 2D-Ising 
model with polynomially more spins. Even more gener- 
ally, one can verify that the partition function of an arbi- 
trary g-state spin model (with finite q), where arbitrary 
pairwise or even /c-body interactions with bounded k are 
allowed, can be written as the overlap between a suitable 
quantum state and a product state. This immediately 
implies that every partition function can be expressed as 
a special instance of the 2D-Ising model. However, in 
general an exponential overhead may be required. 

We further remark that the 2D square lattice does not 
play a special role in this context: there are many other 
models with a similar completeness property [1,3]. For 
example, all Ising models on a graph G whose associated 
graph state \<Pq) is a universal resource for MQC, allows 
one to draw the same conclusions as for the 2D square lat- 
tice. Examples of such other universal models for MQC 
include e.g. hexagonal, triangular and Kagome lattices 
3D lattices as well as 2D lattices with holes. On 
the other hand, all models corresponding to graph states 
\<Pq) which are not universal resources for MQC (in the 
sense of universal state preparation [HI]) are not capa- 
ble of expressing partition functions of e.g. the 2D-Ising 
model (or other complete models). Examples of "non- 
complete" interaction patterns are ID structures such as 
chains or trees, or more generally all graphs where the 
decorated graph G has bounded rank width [ll| . 

7. — Summary. We have established a connection be- 
tween evaluating the partition function of a general class 
of classical spin models and measurement-based quantum 



computation. We have used the universality of the 2D 
cluster states, in particular the possibility of preparing 
any other quantum state by means of projective single- 
qubit measurements from a sufficiently large universal 
state, to show a type of completeness of the classical 
2D Ising model: the partition function of any classical 
spin model (Ising and Potts model on arbitrary graphs, 
and beyond) can be recovered as a special case of the 
Ising model on a sufficiently large 2D square lattice with 
complex couplings. Moreover, we have given an explicit, 
efficient construction of the corresponding 2D model. 

Finally, it is an interesting open problem whether a 
restriction to the real (and thus "physical") parameter 
regime of the 2D-Ising model is possible while keeping 
the completeness property. It would also be interesting 
to investigate how the explicit reductions obtained in this 
paper may be related to previous results regardin g th e 
NP-completeness of the 2D Ising model 0] (see also 12 1). 
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